Parity breaking and scaling behavior in light-matter interaction 
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The light-matter interaction described by Rabi model and Jaynes-Cummings (JC) model is in- 
vestigated by parity breaking as well as the scaling behavior of ground-state population-inversion 
expectation. We show that the parity breaking leads to different scaling behaviors in the two mod- 
els, where the Rabi model demonstrates scaling invariance, but the JC model behaves in cusp-like 
way. Our study helps further understanding rotating-wave approximation and could present more 
subtle physics than any other characteristic parameter for the difference between the two models. 
More importantly, the scaling behavior is observable using currently available techniques, such as a 
superconducting qubit under external driving in circuit quantum electrodynamics. Our idea could 
be straightforwardly applied to the study of Dicke model and spin-boson model. 
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The Rabi model [T[ presents an important prototype 
for the interaction between a single two-level system 
(e.g., a spin) and a quantum bosonic field, which has 
been widely applied to almost every subfield of physics, 
e.g., the cavity quantum electrodynamics (QED) and 
the exciton-photon interaction. Under the rotating-wave 
approximation (RWA), the Rabi model is reduced to 
the Jaynes-Cummings (JC) model Q in which counter- 
rotating effects are neglected by the assumption of large 
detuning and small Rabi frequency. It is generally be- 
lieved that the RWA works worse and worse with the 
increase of the Rabi frequency [3]. So, in the case of 
strong spin-field coupling, the Rabi model, rather than 
the JC model, is required. The differences between the 
two models have been studied by solving the eigenener- 
gies [H, the time evolution Q, the Berry phase Q and 
so on. 

Going beyond the interaction between a single spin-1/2 
and a single quantum mode, the Rabi model has been ex- 
tended to a big-spin (S>l/2) system Q or many spin-1/2 
experiencing a single quantum mode, the latter of which 
is called Dicke model It has been shown that the 
RWA introduced in the Hamiltonian of the Dicke model 
brings about completely different phenomena from the 
non-RWA case in the quantum phase transition || |T(| ■ 
Besides, we may also consider a single spin-1/2 inter- 
acting with a multi-mode quantum bosonic field, called 
spin-boson model [ll], EH > to describe the dissipation of 
a single spin under the bosonic bath. The spin-boson 
models with and without the RWA demonstrate differ- 
ent behaviors of the spin dissipation [13| . 

We focus in the present work on the scaling behav- 
iors in the Rabi and JC models, which could present us 
more subtle physics than the solutions of eigenenergies 



and geometric phases for the difference between the two 
models. The different scaling behaviors are relevant to 
different symmetries and can be understood by parity 
breaking. Specifically, we show the scaling invariance 
only in the Rabi model. In contrast, the scaling behavior 
in JC model behaves much differently with the change 
of some characteristic parameters, such as the detuning. 
This helps us to further understand the RWA that the 
counter-rotating terms significantly influence the scaling 
behavior. More importantly, these scaling behaviors are 
strongly relevant to the dynamics of the spin, which could 
be observed in some experimentally available systems, 
such as the circuit QED system T3- 16|. 

The key point of our investigation is the physics behind 
the scaling behavior, i.e., the parity-breaking induced by 
the local bias field. Our results can be straightforwardly 
used to characterize the Dicke model and the spin-boson 
model although the situations of multi-modes and many 
spins are more complicated. 

We get started from the following Hamiltonian (H = 1) 



(i) 



where A and e are the tunneling and the local bias field, 
respectively, ui and (a) are frequency and the creation 
(annihilation) operator of the single-mode bosonic field, 
and A is the Rabi frequency. a z ^ x are the usual Pauli 
operators for the spin-1/2 and a x = o~+ + <r_ with a± — 
(a x ± ia y )/2. Compared to the standard form of the 
Rabi model, Eq. (1) owns an additional term, i.e., the 
local bias, which does not change the physical essence of 
the model. As discussed later, Eq. (1) connects directly 
to the spin-boson model and to currently experimentally 
achieved systems. 

Eq. (1) can be diagonalized by displaced coherent 
states. The eigenfunction of H s f, has the form 0, [TH, [l9j 
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where c„ and d n are coefficients to be determined later 
and 
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are the displaced coherent states with the displacement 
variable q = A/a;. The complete eigensolution of H s t, 
can be obtained from the Schrodinger equation (See Ap- 
pendix A) [20]. We check the ground-state population- 



inversion (er z ), which is, 



(a z ) = (c -) 2 - (do) 2 



Vn 2 + e-W ■ 



(2) 



where j3 — q 2 , k — e/A, and Cq and <ig~ are defined in 
Appendix A. We first perform the second derivative of 
(a z ) with respect to f3, which yields a reflection point 
/3 C = — ln(2«; 2 )/4, by which Eq. (2) is rewritten as 



(3) 



s /3' 1ii(2k 2 



under the scaling transformation /?' = /3/j3 c . From 
the definition of j3 c , we require k ^= 0, i.e., e =/= in 
the present calculation. For a fixed value of k, the 
population- inversion (cr z ) in Eq. (3) is only relevant to 
/3 , rather than to other characteristic parameters (See 
Fig. 1(a)). So (3 C can be regarded as a scale of the Rabi 
model. In addition, if we set f3 =1, the population- 
inversion (a z ) turns to be a constant — 1/-\/3, implying 
a fixed crossing point with variation of /? . It is more 
interesting to demonstrate the scaling behavior of the 

population- inversion (er z ) with a displaced scaling /? = 



(P-0 ? )/Vvf. Since (a z ((3 )) = -1/ yj 1 + 2e~ 12 ^" , 
which is independent of k under the scaling transforma- 
tion, the population-inversion (a z ) with respect to /3 , 
remains unchanged for different parameters /v, (IS shown 
in Fig. 1(b). 



satisfying 



model [13J , the scaling behavior we show here is only for 
interaction between a single spin and a single mode of the 
quantized field, and is actually resulted from the parity 
breaking regarding the parity operator II = cr x e i7rata . 
If we denote the case of e = in Eq. (1) by H sb , we 
have [H sb ,Tl] = 0, with the ground state of their com- 
mon eigenfunction to be \i/jq) = — ( 

n|^o~) = IV'o")) i- e -j an even parity state of II. The 
even parity breaks down in the variation from e = to 
e 5^ because the Hamiltonian H s b never commutes with 
the parity operator II, i.e., [H s b, II] ^ 0. This parity- 
breaking case plotted in Fig. 1(a) (except the point 
j3' = 0) demonstrates the abrupt variation of (a z ) from 
to -1 at k approaching 0, corresponding to the trans- 
lation from spin non- localization (i.e., superposition of 
the spin- up and spin-down) to spin localization (i.e., the 
spin down). It is a peculiar scaling behavior different 
from those relevant to quantum phase transition in the 
spin-boson model. More interestingly, this sharp scaling 
variation induced by the parity breaking could appear for 
larger k (e.g., k =0.1) under a scaling displacement (Fig. 
1(b)). 
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FIG. 1: (color online) Scaling behavior of the ground-state 
population-inversion, (a) As a function of /3' which has a fixed 
point; (b) As a function of /3" which remains unchanged. 

In contrast to the scaling invariance relevant to the 
critical points of quantum phase transition in spin-boson 



FIG. 2: (color online) (a) The ground-state population- 
inversion for H s b as a function of j3' with k = 10 -6 . (b) 
The ground-state population-inversion for Hj c as a function 
of /?' for A > (the red solid line with cusp- like scaling behav- 
ior) and for negative A (the blue dashed line with no scaling 
behavior) with k = 10 -6 . (c) The cusp- like scaling behavior 
of the ground-state population-inversion for Hj C as a func- 
tion of A/A c with A > for k = 1CT 10 (the black solid line), 
k — 10~ 6 (the blue dashed-dotted line), k — 10~ 2 (the pink 
dotted line), and k = 10 _1 (the red dashed line). For con- 
venience of comparison between (b) and (c), we label A c in 
(b). 



To better understand the parity breaking in the Rabi 
model, we may consider the same treatment for the JC 
model (See Appendix B for the eigensolution and the 
discussion about parity states). In the case of e = 0, the 
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ground-state of the JC model is usually of certain parity, 
and energy level crossing occurs only for A >0. In the 
case of e ^ 0, the parity breaks down. But different from 
the scaling invariance in the Rabi model, the variation 
of (a z ) with respect to ft' behaves with a cusp in the JC 
model (See Fig. 2(a) and 2(b) for comparison between 
the two models). It implies that the spin is always in non- 
localization except the rare cases near the critical point 
A c for localization. For k = 10~ 6 , we have /3 C — 6.73, 
corresponding to A = 2.59lo in the Rabi model and to 
A = 0.09cj in the JC model. These are reasonable values 
for the parameters in the two models. For a more clarified 
demonstration, we plot in Fig. 2(c) (cr z ) versus A/A c , 
where A c = vljA takes values much smaller than /3 C . So 
we could have a zooming-in picture for the variation of 
(<J Z ), which shows strong relevance to k. The cusp-like 
behavior is evidently induced by the parity breaking (i.e., 
for k very close to zero). 

The results above remind us of the role the counter- 
rotating terms playing in the light-matter interaction. 
For the parity operator II with respect to H s (,, we have 
n' = WUU = -cr 2 e i7ra ' a for H JC where U is defined 
in [lTjl- If we neglect the terms regarding e, we have 
[II, H' sb ] = and \U.',H r j c ] = 0. But this does not mean 
that the counter-rotating terms have no extra influence 
on the symmetry in the Rabi model compared to the JC 
model. First, in the case of e = 0, the ground-state of the 
JC model (for A >0) experiences a translation from the 
even parity to the odd in the increase of the coupling A, 
but the ground-state of the Rabi model only stays in the 
even parity. This could be understood in physics from 
whether the energy level crossing happens or not. Sec- 
ond, for e 7^ 0, the counter-rotating (or RWA) effect is 
reflected in different scaling behaviors for spin localiza- 
tion or non-localization in the two models. 

Since (a z ) is observable experimentally, the scaling be- 
havior we investigated above should be an alternative 
way to identifying the RWA effect, which presents more 
subtle and interesting physics. In what follows, we take 
the circuit QED system as an example, in which the 
above interaction can be implemented by a supercon- 
ducting qubit strongly coupled to a microwave resonator 
mode via external driving [16| . The Hamiltonian is given 
in units of ft = 1 by 

Hqed = Y a z+ UJ btfb + G(b<j + +tfa-) 

-r» 2 (e^ t a_+e-^V + ), (4) 

where u q and uj^ are, respectively, the qubit and mi- 
crowave photon frequencies with G the qubit-photon cou- 
pling strength, b (b^) stands for the annihilation (cre- 
ation) operator of the microwave photon. cr Zj ± are usual 
Pauli operators for the superconducting qubit. VLj and u)j 
are the amplitude and frequency of the jth driving field 
(j =1, 2). Assume that the first driving field is strong 
enough, which makes fl± S> ^2, G. In the rotating frame 



first with the driving field wi and then with a large fre- 
quency 3 comparable to fij , we may obtain an effective 
Hamiltonian by setting Q 3 = (uj\ — uj 2 )/2 and neglecting 
fast oscillating terms, 

H eff = ~2 a z H 2 ° X ^ ~ ' 

+ |(6 + 6t)(a + +a_). (5) 

Eq. (5) can be used to simulate the Rabi model and 
JC model by tuning the characteristic parameters. To 
demonstrate the scaling behaviors described above, we 
require f^i — O3 to be exactly tunable from zero to a 
value much smaller than other parameter values. If we 
have fi 2 <C (w& — uji) ~ ^-j ^eff simulates the interaction 
between the superconducting qubit and the microwave 
photon under the Rabi model of Eq. (1), which, in a 
rotating frame, describes Eq. (1), i.e., 

H' ff = Z^<7 g + (Ql ~ n 3) az+{cJb -u; 1 )b%+^(b+bi)* z . 

Alternatively, if O2 ~ w& — uji ^S> § is fulfilled, H e ff turns 
to be the JC model Hamiltian H JC describing interaction 
between the superconducting qubit and the microwave 
photon, equivalent to 

We survey the relevant experimental parameters be- 
low for our purpose. For the Rabi model, to meet the 
condition O2 "C (w;, — uji) ~ G/2, we can adopt follow- 
ing parameters as u q = 2ir x 6.02 GHz, ojf, = 2ir x 6.02 
GHz, oji — 2n x 6 GHz, cj 2 — 2tt x 4 GHz, and a tunable 
G > 2tt x 40 MHz. Besides, we may assume the driving 
fields with the amplitudes O3 = (uji — uj 2 )/2 = 2n x I 
GHz and fJ 2 = 2n x 0.4 MHz, (fii - Sl 3 ) changes from 
zero to 27rx20 kHz, which are realistic values using state- 
of-the-art circuit-QED technology [HI, HI]. While for JC 
model, the result can be obtained by reducing the cou- 
pling strength to G' < 2ir x 4 MHz with other parameters 
unchanged. Thus this circuit QED system can be effec- 
tively tuned to a Rabi model or to a JC model under our 
choice of appropriate values of the characteristic param- 
eters. The observed ground-state population inversion 
(a z ) with respect to (3' or /?" or A/A c defined in Figs. 
1 or 2 is a direct experimental manifestation of physics 
regarding the scaling behavior and the RWA. 

In conclusion, we have investigated the scaling behav- 
iors of ground-state population-inversion in light-matter 
interaction, which is relevant to the parity breaking. The 
different scaling behaviors can be used to identify the dif- 
ference between the Rabi model and the JC model be- 
cause it is more sensitive to the RWA than energy spec- 
trum or geometric phase, and it demonstrates the change 
of the symmetry due to the RWA. We have exemplified 
the circuit QED system to demonstrate experimental fea- 
sibility of our study using currently available technology. 
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The present idea could be straightforwardly extended to 
the study of Dicke model and spin-boson model for pos- 
sible quantum phase transitions induced by the parity 
breaking. 
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APPENDIX 



obvious from the expression of eigenenergies that the 
ground-state energy Eq is smaller than Eq . In the 
case of e = 0, the eigenfunction of the ground-state is 

B. Eigensolution to the JC model 

We consider Hj C = ^fz + %Ox + wa^a + \{aa+ + aV_) 
by neglecting the counter- rotating terms [l7j |. Under 
the condition of a negligible local bias, we may use 
H r j C = ycr 2 + usa^a + \(aa + + aV_) to analyze the 
parity-relevant scaling behavior of the JC model. The 
eigenfunctions of H r j C include |\I>)g — |0)|<?) with energy 
E G = -f and |tf n ) = f n \n)\e) + K\n + l)\g) with 
normalized coefficients /„ and h n to be determined from 



A. Eigensolution to the Rabi model 

Taking |W) into the Schrodinger equation of H s b, 
we have 



[uj(m — q' 
[oj{m — q' 
where D m n is given by 



A 
~2 



2 



= Ed m , 



A x 



D„ 



-2<r 



E 

k=0 



(-1) 



_ k VmM(2q) m+n - 2k 
(m - k)\(n - k)\kV 



In our case with the condition A/cj <C 1, the 
terms of D m ^ n with m ^ n play negligible roles in 
the equations compared to other terms with m = n. 
So the equations above can be simplified to the case 
remaining the terms of m = n, which yields follow- 
ing analytical solutions, that is, the eigenenergies 

E± = uj(m - q 2 ) ± Je 2 + A 2 D 2 a ~/2, and the coeffi- 



cients c± = Mm/v 1 + (Mm) 2 and d± = l/ x /l + (^m) 2 
with /i± 



e ± 



A 2 D 2 ' 

^ m , rn , 



/(AD r . 



It is 



(nco + —)f n + Xy/n + lhn = Ef n , 



[(n + l)w - —]h n + XVn + If a = Eh r , 



We may obtain E n> ± 
\ v/(w- A) 2 +4A 2 (n+l) with 



w(n + |) ± 



|*)„,+ = cos -f\n)\e) + sin f\n + l)\g), 
|*) n _ = sin y|n)|e) - cos y|n + l)\g), 



and cos( 



A- 



As we focus on the 



^(^-A)2+4A2(„+l)' 

ground state, the possible lowest state should be | o. 

with the energy E ,- = [w- - A) 2 + 4A 2 ]/2. 

If E ,_ = E G , i.e., [u - ^(uj - A) 2 + 4A 2 ]/2 = -A/2, 
we have the level crossing between \^)g and | x E')o,- at 
A c = \JuiA with A > 0. For A/A c < 1, the ground state 
is Y$)g — 1 0) | e?) of even parity because n'|\&)(j = \^)g- 
In contrast, if A/A c > 1, the ground state is | v I')o,- = 
sin^|0)|e)-cos^|l)|5) of odd parity due to n'|*) ,_ = 

— | SI/) o, - The parity breaks down if e is tuned from zero 

to non-zero. 
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